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[x^ . Abstract 



The complete soft-enhanced and virtual-gluon contributions are derived for the quark coefficient 
functions in semi-inclusive e+e annihilation to the third order in massless perturbative QCD. 
These terms enable us to extend the soft-gluon resummation for the fragmentation functions by 
two orders to the next-to-next-to-next-to-leading logarithmic (N^LL) accuracy. The resummation 
exponent is found to be the same as for the structure functions in inclusive deep-inelastic scattering. 
This finding, together with known results on the higher-order quark form factor, facilitates the 
determination of all soft and virtual contributions of the fourth-order difference of the coefficient 
functions for these two processes. Unlike the previous (N^LL) order in the exponentiation, the 
numerical effect of the N^LL contributions turns out to be negligible at LEP energies. 



Semi-inclusive e+e^ annihilation (SIA) via a virtual photon or Z-boson, e^e^ ^ y/Z ^ h +X, 
is a classic process probing Quantum Chromodynamics (QCD), the theory of the strong interaction. 
A wealth of precise measurements have been performed, at various center-of-mass (CM) energies 
of the total fragmentation function 

where h stands for a specific hadron species or the sum over all (charged) light hadrons, see Ref. [1] 
for a general overview. In the CM frame the scaling variable x is the fraction of the beam energy 
carried by the hadron h, and = s h the square of the four-momentum q of the intermediate 
gauge boson. In perturbative QCD, the total (angle-integrated) fragmentation function F/^ = F^, 
as well as the transverse {Fj), longitudinal {Fi) and asymmetric {Fa) fragmentation functions for 
the double-differential cross section do'^ /dxd cosQh [2], are given by 

Here are the parton fragmentation functions, the final-state (timelike, = q^) analogue of 
the initial-state (spacelike, = — <?^) parton distribution functions in deep-inelastic scattering 
(DIS). Without loss of information in the present context, the renormalization scale of tts and the 
factorization scale of have been identified with the physical hard scale in Eq. (O. The 
coefficient functions ^ are defined via expansions in the strong coupling = a^/ (47r). 

Here we are interested in the dominant (anti-) quark contributions to F/', Fj and F^, 

C^q{x,a^) = Oew(5(l-Jc) + asci^W + ^^s^ci^W + a^c^^Jiix) + ...). (3) 

The electroweak prefactors Oew can be found in Ref. [2]. The first- and second-order coefficient 
functions have been calculated long ago in Refs. [3] and [4], respectively. More recently the latter 
results have been confirmed (and some typos corrected) in two independent ways in Refs. [5,6]. 

(3) 

The three-loop corrections Ca {x) have not been derived so far. 

The coefficient functions in Eq. ([3]) include large-x (threshold) double-logarithmic enhance- 
ments of the form a" ( 1 —x)^^ln^{l —x) with = 0, . . . , 2n — 1 . Such contributions, which spoil 
the convergence of the perturbation series at sufficiently large values of x, can be resummed by the 
soft-gluon exponentiation [7, 8]. For the process at hand this resummation has been worked out to 
next-to-leading logarithmic (NLL) accuracy in Ref. [9]. The inclusion of this resummation has led 
to improvements in a recent global fit of fragmentation functions [10]. Hence an extension of the 
soft-gluon exponentiation for e+e^ ^ y/Z ^ h +X to a higher accuracy is not only of theoretical 
but also of phenomenological interest. 

In this letter we employ the analytic continuation approach of Ref. [5] to derive the soft and 
virtual contributions to the third-order coefficient functions in Eq. (|3]). These results are then 
used to extend the results of Ref. [9] to the next-to-next-to-next-to-leading logarithmic (N^LL) 
accuracy reached before for inclusive deep-inelastic scattering [1 1] and the total cross sections for 
lepton-pair and Higgs-boson production in proton-(anti-)proton collisions [12, 13]. A substantial 
intermediate step towards the present extension has been taken before in Ref. [14]. 
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Up to small contributions from higher-order group invariants entering at the third and higher 
orders, the soft plus virtual contributions are the same for the DIS quark coefficient functions 
for Fi, F2 and F3 [15, 16]. The same holds for the corresponding (in this order) SIA coefficient 
functions for Ft, Fi and Fa- Hence we will drop the index a from now on, and refer to the former 
coefficient functions collectively as c)^ {x), and the latter as Cj\x). 

In this limit the bare (unrenormalized and unfactorized) partonic DIS (spacelike) structure 
function F^ is given by [17, 18] 

F^{alQ^) = 5(1 -x) + l^ia^y (^) F,^^, (4) 

z=l VA* / 

with 

F^, = 2iFi5(l-x)+5i 

F^^ = (2j2 + (Ji)')5(l-x) + 2iFi5i+52 

F^s = (2j3 + 2iFiiF2)5(l-x) + (2y2 + (Jl)')5l+2Ji52 + 53 • (5) 

Here /j is the scale of dimensional regularization with D = 4 — 2e, and the bare strong cou- 
pling. J/ represents the Z-loop quark form factor [17-22]. The .x;-dependence of the real-emission 
functions Sk is given by the D-dimensional -i-distributions 

AeW = [{l-x)-'-'']+ = -l5(l-x)-f£^^©, (6) 

where we have introduced the abbreviation = [(l~-'c)^ln^(l— jc)] + . 

The transition to the bare SIA (timelike) fragmentation functions F^ is performed as follows: 
In Eq. ^ the factors 2^/ are replaced everywhere by 2ReJ/'^ and all products J^t^/ by 
where is the complex /-loop timelike form factor which can be obtained from the spacelike [Fi 
by Eq. (3.3) of Ref. [17]. The analytic continuation of the real-emission terms Sk is carried out as 
discussed in Ref. [5]. In fact, these functions turn out to be the same for the spacelike and timelike 
cases (this holds only in the present large-x limit, not for the full real emission contributions). 
Finally the standard renormalization and mass factorization is performed to the third order for the 
resulting timelike analogue of Eq. dS]), yielding the 'D^ and 6(1— jc) terms of Ca (x) in Eq. ([3]). 

For the convenience of the reader, we include also the large-x limits of the well-known first- 
and second-order MS coefficient functions [3,4]. As expected from the above discussion, these 
and the third-order coefficient function share all non-^2 terms with their spacelike counterparts, 
hence we will present them via the corresponding differences 5-pg c„ = Cj^ — c^^\ The results read 

djsc.ix) = 12^2C^6(l-x) , (7) 

Sts^zW = 48C2C2©i-36C2C|2)o 

+ { (- 108 + 24 ^2) C2 + ( ^ - 24 ^2 ) QC^ - y C^n^ } ^2 6( 1 -x) , (8) 



2 



- { (324 + 96 ^2) C3 - ( ^ - 192 ^2 ) C^Cj + ^ Cjuf } ^2 
+ { (306 + 2 16 ^2 - 96 C3) - ( - 248 C2 - 480 C3) QC^ 

l^-32C2)c2n^}C2 2^0 

993 / 13457 220 



+ 



+ {( i!^ + ISOfc -936C3 + 72a)c3 -( lili^ + ^ - 1616^3 
108,,\_„2 /74728 528 _2 



+ ( f + f C2 - 80 C3) C|„, - ( + ^ C2 - 96 C3)qC,„, 

+ (^ + yC2)Q«;}C28(l-x). (9) 

Here Ca and Cf are the standard group invariants, with Ca = 3 and Cf =4/3 in QCD, and 
the number of light flavours. denotes Riemann's ^-function. The third-order SIA coefficient 
functions can be obtained by adding the corresponding DIS results given in Eqs. (4.14) - (4.19) 
and Appendix B of Ref. [15], see also Eq. (3.8) of Ref. [16]. The first half of Eq. ^ agrees with 
the result of Ref. [14], the 5(1 —x) contribution in the second half has not been presented before. 

Below we will need the A'^-independent parts = 5ts c^{N) of the Mellin transforms 

of Eq. © - dll) obtained via 

= j^dx [x^^^ - 1) a{x)+ (10) 
together with 5(1 —jc) 1. These contributions are given by (Ye is the Euler-Mascheroni constant) 

(11) 

' -C2(l08-48C2-36Ye-24Y2)- — Qn^, (12) 



^2 


^ Sts.?oi 


= 12Cf 


^2 


^ ^TSS02 


= Cj^Cf 


^2 


^ ^TS.?03 





(— + 18 C2 - 792C3 + — C2 - 306Te + 288TeC3 - 162j^ 
+ 96YK2 + 12J^ + 24 J^) + C,C| ( - Hl^ + 482^2 + ^ & 
- + ^Te - 160TeC2 -480T.C3 + ^Te^ -96y?C2 

+ ^,|)^CiC,(Zf?-,96C2-.056C3 + f C?) 

^2 /667 272 ^ 1672 ^ 268 2 16 3 

+ C|n^(^ — -44^2- — C3--^Ye + 16YeC2-— Ye -yYe 

/ 23504 16 ^ ^^^\ ^ 2/ 1624 16 ^ \ 
+ CAC^nf(^~-^ _^2 + 96C3)+C^n2(^^ + -^2) • (13) 

The corresponding DIS coefficients can be found in Eqs. (4.6) - (4.8) of Ref. [11]. 
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For processes such as DIS and SIA, the dominant large-x/large-A^ contributions to the MS 
coefficient functions can be resummed by a single exponential in Mellin space [7] 

C^(e') = ^o(e') -expiG^le^)] + o{N-Hr^^N) . (14) 

The prefactor collects, order by order in the strong coupling constant tts, all A^-independent 
contributions. The exponent contains terms of the form In^A'^ to all orders in ttg. Besides the 
physical hard scale 2^ (= =f^^ in DIS/SIA, with q the four-momentum of the exchanged gauge 
boson), both functions depend on the renormalization scale /j^ and the mass -factorization scale jU^. 

The exponential in Eq. (fT4l) is build up from universal radiative factors for each initial- and 
final-state parton p, Ap and 7p, together with a process-dependent contribution A'"^ The resumma- 
tion exponents for DIS and SIA [9] take the very similar form 



'DIS 

-A' 
'SIA 



InAq + ln7q + InAj^Jg , 
InAq + ln7q + InA^'^ • 



(15) 



The radiation factors are given by integrals over functions of the running coupling. Specifically, 
the effects of coUinear soft-gluon radiation off an initial-state or 'observed' final-state quark are 
collected by 



lnAq(e2,A'/) 



r- '-1 r('-^)'Q'dq' 







1 .,A'-1 

dz ^ 



(16) 



CoUinear emissions from an 'unobserved' final-state quark lead to the so-called jet function. 



ln7q(2^ 



dz- 



1 



1 



dq^ 



A{a,{q^))+B{a,{[l-z\Q^)) 



'(i-z)2e2 r 

Finally the process-dependent contributions from large-angle soft gluons are resummed by 

1 



(17) 



InA^^Xe 



1 ^^-1 



dz- 



■D{a,{[l-zfQ^)). 



1-z 



The functions in Eq. (fT4l) and A, B and D in Eqs. (fT6l ) - (fTSl ) are given by the expansions 



at 



l=lQ l = lo 



(18) 



(19) 



where lo = with = 1 for F = ^q, and Iq = \ else. 

The known expansion coefficients of the cusp anomalous dimension (the coefficients of 
1/(1 —x)_^_ in the MS quark-quark splitting functions) read [23, 24] 

Ai = 4Cf 
A2 = SCf 



16C, 



C, 



( 



245 



V 24 
+ CArif - 



67 „ 11 ^ 11 ;-7\ ^ / 55 ^ ^ \ 



209 10 ^ 7 . 

\ C2 C3 

108 9 ^ 3 ^ 



+ nti 

27 



(20) 
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The first three coefficients of the jet function (flTl) are given by [7, 1 1, 25] 



Bi 
Bi 



-3C 



F ■> 



ct 



+ Cpn 



2 
247 

~7n 



+ 12^2-24^3 



3155 44^ ■ 



(21) 



(22) 



53 



Ci 



29 

y-18C2 



68^3 



288 



+ 32 C2C3 + 240^5 



712 

■46 + 287 ^2--^C3 



272 



+ ClCp 



+ Cpnf 



599375 



3 

32126 



Cl- 16^2^3-120^5 



729 



^ 21032 ^ 652 ^2 

+ -T^C2 + -;:^C3--r^C2 



176 



81 



5501 



+ Cj^CpH 



54 
160906 



50C2 + yC3 



729 



9920 



^2 



27 

+ Cpnj 
776 



15 



8714 232 ^ 

1 C2 ■ 

729 27 



^2^3-232^5 
32 . 



27 



r 208 ^ 

^3 + ^^: 



(23) 



Together with Eqs. (fTTI) - (fT3l) . all functions but D in Eqs. (fT6l) - (fT8l) are known to order a^^. 
Consequently the first three coefficients of D^^^ can by determined by comparing the ttg -expansion 
of Eq. (fT4l) with the fixed-order results (|7]) - This procedure yields 



SIA 







(24) 



for = 1, 2, 3, hence A™/^ 



1 to at least N^LL accuracy. Di—0 was, of course, included in the 
NLL resummation of Ref . [9] . However, B2 was unknown at that time, and only B2 + £>2 could be 
extracted from the two-loop results of Refs. [4] alone. 

As expected from the identity of the DIS and SIA soft-emission functions Sk in Eq. ([5]), there 
is a strong similarity between the respective coefficient functions also in the framework of the 
soft-gluon exponentiation - recall that 



D 



DIS 



0, 



A int 
^DIS 



1 



(25) 



was proven to all orders in tts in Refs. [26, 27]. We expect that such a proof can also be derived 
for SIA. For the time being assuming the all-order validity of Eq. ((24)) . the difference between the 
SIA (timelike, T) and DIS (spacelike, S) large-A^ coefficient functions exponentiates as 



5TsC^(e') = 5Ts^o(e')-exp[G^(e^)] + o(^-Mn"A^) (26) 
where, after performing the integrations in Eqs. ([T6l ) - ([T8] ). the function takes the form 

G^(e2) = \nN-g,{X)+g2{'k)+a,g^(X)+aUAi^)--- (27) 

with X = Po In A^. The first three expansion coefficients of 5jg for = = Q have been 
given above in Eqs. (fTTI) - (fT3]) . We will address the fourth-order coefficient below. 
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The functions to have been derived in Refs. [7, 1 1, 28, 29]. For completeness we include 
these functions, also here restricting ourselves to choice /j^ = /j^ = 2 of the scales: 

^piS(^) _ Ai(l-ln(l-X)+X-Un(l-X)) , (28) 

gfS(X) = (Aipi-A2)(X + ln(l-X)) + ^Ai(3iln2(l-^) 

-(Aiye-5i)ln(l-X) , (29) 

8?'H^) = ^(Ai|32-Aip2+A2|3i-A3)(l+X-^) 

+ (AiPiYe +A2P1 -5i|3i) (1 - - ^4^) 

-(Aip2 + ^Ai(Y2 + C2)+A2Ye-5iye-52)(l-Y^) , (30) 



and 



In(l-X) 



-Aipi (Ye" + ^2) +A2pf - 2A2piYe - A3P1 + 25i|3iYe + 252pi) , . , , 2 

(1 -A)^ 

- (AiPl -AiPi|32)^^3^ + (^Aip3 _A,M2 + ^AiPs) ln(l - A) + (AiPf 
-Aipip2-Aip?Ye+Ai|32Ye-A2p?+A2p2 + 5ipl-5ip2)(^-Y^ 

+ ^ (T^) + \ il^'^' " ^'^^'^^•^^ ~ l^'^' " ^^l(3YeC2 + Ye +2^3) 

+ A2plYe - A2(Ye + ^2) " ^^2(3? + IA2P2 + ^Agpi - AgYe " -52|3i 

D i D i 

+ 5l(Ye+C2)+252Ye+53)(l-— ^)+^(Ai|33-2Ai(3ip2+Ai|33 



(1-A)^ 

+A2p2-A2p?+A3pi-A4)A. (31) 

Factors of Po = 11/3Ca — 2/3ny have been suppressed in Eqs. (l28l) - (|3TI) for brevity. The 

dependence on Po is recovered by A^/^q, Bk Bk/^Q, % ^k/^o^^ and multiplication of 
g2 and g4 by Po and Pq , respectively. Note that Eq. (|3TI) includes all known coefficient of the beta 
function of QCD, see Ref. [30] and references therein. 

All parameters entering Eqs. (1281) - (|3T]) are known except for the four-loop cusp anomalous 
dimension A4. The small (see below) impact of this quantity - which first occurs in the In-^A'^ 
contribution to 6jgC^ - can be included by a Pade estimate as in Ref. [11], backed up by a 
recent calculation of one Mellin moment of the fourth-order quark-quark splitting function [31], 
cf. also Ref. [32]. E.g., for n^ = 5 one may use A4 ^ 1550 (recall our small expansion parameter 
as = tts/ (47r)) and assign a conservative uncertainly of 50% to this value. 



Due to the vanishing of 6^5 ^qq the two highest logarithms, a^ln^^ N and al.\n^^^^ N, are the 
same for the SIA and DIS structure functions to all orders in tts. The expansion of Eq. (|26l) 
with Eqs. (fTT)) - (fT3l) provides the six highest logarithms, cf. Ref. [11], of the coefficient-function 
difference 5-pg C^, ai. In^^^^N with a = 2, . . . , 7, at all orders from the fourth. In particular, all InN 
enhanced terms are thus fixed at order a^. After transformation to x-space these contributions read 

5tsC4W = 96C2C>5-{360C4 + ^QC3-i^c3n^}^2©4 



3 3 

- I (432 + 576 C2) - (3552 - 576 C2 ) QC^ + 576 Cjn^ - ^ Cfc^ 
+ ^QC2nf-^c2n2|^2 2)3 + {(1674 + 2160C2 + 192^3)C4 
_ ^ 25238 _ _ 2330 j ^^^3 + 1^ _ 352 C^n^ 

-(^-528C2)cic,^+(^-96C2)QC,V-^C,V}C2^2 
+ {(ll22 + 936C2-4320C3-^ClK^-(^ + ^C. 

- 3584^3 - ^ C|)qc3 + ( If + If ^2 + 64C3)C|., 

/ 224230 17176 ^ , 5184 . . 

+ ( — ^ ^ C2 - 7392 C 3 + ^ Cl) CicI 

/ 69728 3056^ \^ ^2 /4888 64 ^ \ 2 2U 

- ( —9 ^ ^2 - 576 ^3 j QC + _____ C 1 ^2 © 1 

- { ( ^ + 33 1 2 ^2 - 3288 ^3 + 792 C 2 + 1 92 - 5 1 84 ^5) 
-(^ + ^^2-8816^3- 1452^1-1728^2^3 

-1440C3)QC3^(^-^C2-l^C3-.^Cf 

- . 1- f . t- \ ^7^9 / 3551 13568 „ 688 \ 

+ 1664^2^3 +2784^5 C^Cl + -— + — — C2 + ^3 



^^-r^-S2^ — S3;-F'y 

^ 1983208 66392 .oo.. ^ ^^^^ 2^ 
- ( ^ ^2 - 2336 ^3 + ^ ^2 j CACpHf 

( 135020 464^ 128 ^ \ 2 2I ^ 
+ (^^^ ^C2 + ^C3)c|n2|^2 2)o + ... . (32) 

The first four terms correspond to a NNLO + NLL accuracy as first obtained for DIS in Ref. [33]. 
For the present case these terms have been presented, in a different notation, already in Ref. [14]. 
The coefficients of and (recall the definition below Eq. ([61)) are new results of the present 
study. The latter coefficient depends on our assumption that Eq. (|24l) extends to = 4. 

The fourth-order result (|32l) can be verified, and extended to the 5( 1 —x) contribution, in the 
following manner. Eq. ([5]) is extended to the fourth order, 
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S,4 



(254 + 2 J3 + (^2) ) 5( 1 - + ( 2iF3 + 2 iF2) 5 1 



+ (252 + (^l) )52 + 2ji53+54 



(33) 



and is subtracted from its timelike counterpart obtained as discussed above. Assuming that also 54 
is identical in the two cases, the only unknown in ^j^^A order is the four-loop anomalous 
dimension A4. All other unknown quantities, such as the and contributions to the spacelike 
three-loop form factor [18, 21, 22] (also the latter new result is not needed in the present context), 
drop out in this difference. Also the four-loop form factor is known from its exponentiation [34] to 
a sufficient accuracy in £ [18]. The soft and virtual contributions to 5jgC4 are then extracted from 
the fourth-order mass factorization formula (here given in terms of the bare coupling) 



5js^4 



StsC4 + 



^PoPi-^PiPo-^PoPi + ipo^i 



+ 



+ 
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PoPo'+PoPo-r^d 







Sts^i + 





6xs«2 



3po2-^/'oPo + ^Po^] 5ts^2 + [3Po-/'o]5ts«3 + e-terms . 



(34) 



For brevity we have suppressed the . . . terms which form a consistency check but do 
not provide new information. The functions a„, b„ and d„ are the £\ and £^ contributions, 
respectively, to the D-dimensional coefficient functions at order a", cf. Ref. [16], and denotes 
the N"LO quark-quark splitting functions. In jc-space obviously all products of these functions in 
Eq. (|34l) have to be read as Mellin-convolutions. 

The determination of 5 C4 from Eqs. (1331) and (l34l) reproduces the result in Eq. (|32l) — hence 
D^iA ^ ^Dis -j^ gq corresponds to dj^Sk = in Eqs. ([33]) and their higher-order 

generalizations 

^2 ^ ^TS Q 



5(l-x) 



and includes the final large-^: coefficient, 
'''' 3779C,-3816C3- 



+ 



+ ■ 



2 

3 1856 
105~ 
62452 



1216a)C; + 



5 

191411 



+ 4080^2^3 + 14880^5 
153802 



12 



+ 



+ 



9 

14817221 
324 



^2 8128 ^ ^ 67328 ^ 

C2 + -^C2C3-— 



9 

102472 
105 



C2- 42808^3 
C| + 4064 ^|)cjQ 



63347 ^ 1856680^ 5306 ^2 



+ 6256^5 + ^^1 



21 



992^^ 



27 



45 



2032^2^3 



13294462 206162 , 



243 



416032 



^3 - 1 100 + 1 936 C2C3 + 8976 ^5 C^Cj 



409 23350 



^2 + 6840^3 



/ 706405 187834 ^ 



55592 
45 



27 



2272 ^ ^ 6272 

-^^2^3 + -^ 



416384 ^^^ 6932 ^2^3^^^^^^ 



27 



45 



8 




(35) 



See Ref. [15] for the flu diagram class leading to the term with d dabc / = 5/18^^^ in QCD. 
The numerical effect of this contribution is very small and will be disregarded in the following. 



hence (up to the residual uncertainty due to A4) the seventh tower of large-jc logarithms from order 
for this difference. For nj: = 5 quark flavours, the numerical expansion of 5jg^Q is given by 

Sts.?o(«s) ^ 2.094as(l + 1.463as + 2.749as^ + {6.659 + 0.094A4/1000}as^ + ...) . (36) 

Thus the two new terms form a correction of almost 5% at tts = 0. 12, with a negligible uncertainty 
from the missing exact value of A4, and the fourth-order contribution is less than half of the previ- 
ous term for ttg < 0.2. It is well-known that the coefficients in Eq. (|36l) are due to ^2-terms (i.e., 
powers of Ti^) from the analytic continuation of the form factor which are subject to a separate 
exponentiation (see, e.g., Refs. [34]). The corresponding results for the SIA and DIS cases read 



The pattern of the corrections in Eq. (|37l ) and the size of the a^-term in Eq. (|36l ) strongly suggests 
that the fourth-order contribution to q amounts to less than 0.5% for ttg = 0. 12. 

The coefficients of the known In^A^ terms are given in Table 1 to the tenth order in tts, using 
the notation c^^ for the coefficient of a^Xn^^^^^^ N in C|[^. Hence, as in Ref. [11] for the DIS 
case, the coefficients of the leading (next-to-leading etc) logarithms are denoted by c^j (c^j ^^^)- 
The qualitative pattern of these coefficients is similar to the DIS case (where all numbers 
are smaller). The higher-order coefficients rise very rapidly, by about an order of magnitude or 
more, with a until a = k — Bk4 without showing the larger-a turnover of the DIS coefficients, cf. 
Table 1 of Ref. [11]. Indeed, the coefficient for the two cases are very similar for a <^ k, but the 
SIA coefficient are more than double their DIS counterparts ata>k where the numbers are large. 

Consequently the higher-order soft plus virtual contributions are qualitatively similar, but larger 
in the timelike case. The numerical size of its resummed coefficient function ([141) is illustrated in 
Fig. 1 for a value of tts corresponding to LEPl, s = MJ. Obviously the size of the coefficient 
function, as well as the relative impact of the new N^LL and N^LL corrections, increases towards 
lower CM energies. Nevertheless one can conclude from Fig. 1 that the accuracy now reached for 
the dominant large-x/large-A/^ contributions should be sufficient for the foreseeable future. 



The Mellin transform of these equations provides the prefactor 5jg^Q4 in Eq. (|26l) . and 



^T,o(«s) 
^S,o(«s) 



1 + 1.045as + 2.266 + 4.703 a^^ + ... , 
1 - 1.050as - 0.797 - 1.056as^ + ... . 



(37) 
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K 




^k2 


^k3 


^k4 


^k5 


^kb 


n /in 


1 


2.66667 


7.0785 


— 


— 


— 


— 


— 


2 


3.55556 


25.6908 


105.621 


104.34 


— 


— 


— 


3 


3.16049 


43.3408 


309.335 


1016.50 


2306.0 


2090 


— 


4 


2.10700 


46.6020 


514.068 


3125.96 


11774.1 


23741 


4664 


5 


1.12373 


36.4525 


577.143 


5393.82 


32365.2 


110255 


29009 


6 


0.49944 


22.3131 


481.110 


6314.54 


55037.7 


293931 


119399 


7 


0.19026 


11.1933 


315.972 


5515.83 


65426.2 


506294 


294105 


o 

8 


0.06342 


4.7503 


170.251 


3808.07 


58765.0 


1 C^C\ AC\ 

618949 


487117 


9 


0.01879 


1.7455 


77.500 


2160.26 


41980.1 


574684 


589591 


10 


0.00501 


0.5652 


30.470 


1035.7 


24725.4 


425171 


551698 



Table 1: Numerical values of the five-flavour coefficients c^^ of the ln^^^"+^A^ contributions 
to the coefficient function C^^. The first six columns are exact up to the numerical truncation, 
and the same for iv, Fj and Fa- The seventh column neglects the tiny (and non-universal) fl\\ 
contributions, and uses the estimate A4 = 1550 for the four-loop cusp anomalous dimension. 
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Figure 1 : Left: the LL, NLL, N^LL and N^LL results for the threshold resummation (fT4l) of the SIA 
coefficient functions Q in A^^-space. Terms to order a" are included in q for the N"LL curves. 
Right: the convolutions of these results with a schematic large-jc shape for the quark fragmentation 
functions, using the standard 'minimal prescription' contour [8] for the Mellin inversion. 
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To summarize, we have first employed the close relation between the perturbative corrections 
to the structure functions in deep-inelastic scattering (DIS) and the fragmentation functions in 
semi-inclusive e+e^ annihilation (SIA), see also Refs. [35], to derive the complete soft and virtual 
corrections to the third-order quark coefficient functions for the latter observables. 

This result then made it possible to extend the soft-gluon exponentiation in SIA from the next- 
to-leading logarithmic (NLL) contributions [9] by two orders to N-^LL accuracy (we confirm the 
intermediate results in Ref. [14]). It turns out that the resummation exponents are the same, pre- 
sumably to all orders, for the DIS and SIA coefficient functions. Hence the threshold enhancement 
is structurally identical in the two cases, and the same thus holds for the class of large-x 1/2^ 
power corrections associated with the renormalon ambiguity of its perturbation series [27, 36]. 

The N^LL exponentiation fixes the seven highest large-jc logarithms at the fourth and all higher 
orders in ttg. The especially simple connection between the soft and virtual contributions to the DIS 
and SIA coefficient functions also facilitates a full N-^LL resummation of the SIA — DIS difference, 
including the next-to-next-to-next-to-leading order a^5(l — contribution to this difference. 

Since the prefactor of the resummation exponential is larger in SIA than in DIS, the soft-gluon 
enhancement is numerically larger in the former case. However, while the N^LL contributions are 
still significant at LEP energies, the N-^LL corrections are practically negligible, indicating that a 
sufficient perturbative accuracy in the large-x limit has been reached with the present results. 
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